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, ^ ■ We calculate the volume corrections to the two-pion matrix element of twist two operators. 
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1. Generalized parton distributions (GPDs) provide a unified parametrization for nonforward 
fiadronic matrix elements [H [2], [Sj H] . Understanding GPDs in detail is therefore tantamount 
to understanding in large parts the internal structure of hadrons. The comprehensive review 
of experimental programs as well as theoretical works can be found in Refs. [5], El [7] . 

As a rule GPDs enter into physical observables only within convolutions. These are rel- 
atively simple at leading order in «<, but become increasingly complex at higher orders, see 
e.g. [8]. It makes the extraction of GPDs from experiment a highly non-trivial task. There are 
much hopes that the lattice QCD can be helpful for solving this problem. There already exists 
a number of lattice results for moments of GPDs, see [9l [lOl [H] and references therein. 

The lattice evaluation of the moments of GPDs is based on calculation of the form factors 
of local operators, very similar to the case of the usual electromagnetic form factors ^12] . 
Unfortunately, at present the lattice calculations with dynamical quarks can be done only for 
unphysically large quark masses and thus pion. The mass of the pion affects strongly the 
spatial extent of the nucleon and hence its form factors. Since the calculations are restricted 
to finite size lattice additional corresponding corrections have also to be taken into account. 

Both problems - calculation of the mass and volume corrections - can be addressed in 
the framework of chiral perturbation theory (ChPT) fl3[ [HI |15l [16j. These corrections were 
already calculated for a number of physical observables (masses, decay constants, etc.). In 
particular, the mass corrections to GPDs (pion and nucleon) have been calculated recently 
in Refs. fTfi [T8| [T9l [201 [2T] . The calculation of the volume corrections in the framework of 
ChPT was initiated by the paper of Gasser and Leutwyler [16] . Recently, the problem became 
rather actual due to the needs of lattice QCD. A review of the present state of the art and 
references can be found in Ref. [23] (see also Ref. [22| where the pion GPDs were calculated in 
partially-quenched ChPT at finite volume). 

On a lattice the overall spatial symmetry is reduced from the usual Lorentz group to the 
hypercubic group. As a result the restrictions imposed by this symmetry on the form of matrix 
elements of twist -two operators are less stringent than in the continuum limit. The number 
of independent formfactors increases very fast with the spin of an operator, and the reduced 
symmetry also results in additional mixing among operators [23|. All this makes the extraction 
of the GPD moments from lattice data a rather complicated problem. Therefore, the knowledge 
of volume corrections to the matrix elements of twist-two operators would be very helpful for 
determining GPDs moments on the lattice. In the present Letter we report the calculation of 
the volume corrections to the lowest moments of vector and tensor GPDs for the pion. 

2. To set the notations we recall the main ingredients of chiral perturbation theory [T3l [Ti] . 
ChPT is an effective theory describing pion interactions in the limit g, -C A^, where g is a 
generic momentum and A^ is the scale of chiral symmetry breaking. We will use the standard 
(9(g") power-counting of ChPT. The Goldstone boson fields are collected in the matrix-valued 
field U. The leading-order pion Lagrangian reads [Tl| 

4^2 = ^ tr {D.UD^V^ + x^U + xU^) , (1) 
where the covariant derivative and the external field tensor x are defined by 

D^U = d^U -ir^U + iUl^, X = 2Bois + tp) . (2) 
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Here s and q denote external scalar and pseudoscalar fields, which count as quantities of 
order O(g^). Furthermore, and are external right- and left-handed vector fields with 
intrinsic chiral power 0{q). The two leading-order parameters appearing in ([1]) are the pion 
decay constant F (normalized to F 92 MeV) and the two-flavor chiral condensate Bq = 
— (0|gg|0)/F^, both evaluated in the chiral limit [T?]. Throughout this work we use the non- 
linear representation for the pion fields vr", 

U = exp{iT^7ry{{7r/Ff)/F}, (3) 

where r"^ are the Pauli matrices. The choice of the function ip is more or less arbitrary. We 
find it convenient to fix it by the condition dfi{U) = Y[t=i dr^"", where dfi{U) is the invariant 
measure on the group SU{2). This requirement results in the following equation 

^{x'') - ^ sin (2x ifix^)) = , ^{x^) = 1 + + Aa;4 ^ ^(^6) (4) 

The leading order pion Lagrangian in the parametrization ([3]), (jlj) reads 

= \ {{d^ - ^'-'} + {H-dnf - Adnr} - ^^{nr + O(vr^) . (5) 

The specific feature of this parametrization is the absence of pion field renormalization at one 
loop, = 1. 

The chiral perturbation theory for pions in a finite volume was formulated in Ref. [T6]. It 
was argued that imposing periodic boundary conditions on the fields, 

'7T°'{Xo, Xi + L, X2, X3) = 7r"(xo, Xi, X2 + L, X3) = 7r"(xo, Xi, X2, X3 + L) = 7r"(xo, Xi), 

does not affect the Effective Lagrangian ([T]) 0. The only effect is a modification of the pion 
propagator according to 

G"^'(x, x') = J2 Gti^o -x'o,x-x' + nL) , (6) 

n 

where the sum is taken over integer n and Goo is the pion propagator in infinite volume 

t J [ZTTpp'^ — m^ + tu 

The imposed boundary conditions reduce the Lorentz symmetry to Rq x H{3), where Rq is the 
reflection xq — » — Xq and H{3) is the cubic group. 

3. The first moment of the isovector GPD for a pion, (see Ref. [18] for a definition) is related 
to the matrix element of the isovector vector operator = qT""y^q 

{7i\p')\Q';,\7r'^{p)) = ze'^'^Ai=\t)P,, Ai=\t) = ^dxH'=\x,i,t) . (8) 

^We will not require a periodicity in time direction. 
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Figure 1: One- loop graphs contributing to two-pion matrix elements of the pion operator. The 
operator insertion is denoted by a black blob. 



Throughout this Letter we use the standard notations for the kinematical variables 
P = lip + p'), A = p'-p, t = A\ ^ = -^' 

where A+(P+) is the light-cone projection of the vector A(P). 

In ChPT the operator takes the following form (we will follow the notations of Ref. |18j ) 

0; = -'-F'{L'; + R;)+0{q'), (10) 

where = Wd^U, = Ud^W and = L'^^t'^ and the same for R^. Calculating the 
corrections to the matrix element of this operator in ChPT one determines the formfactor 
A{"^{t). The infinite volume result for A{"^{t) at order O(g^) reads [H] 

A,^,{t) = 2(l + Fy{t)-ll{^i)t/F'), (11) 

where is a low energy constant (LEC), 

J(t) = 2 + crlog^, (12) 
cr + 1 

= AttF and cr = a/1 — Am? /t. To obtain this result one has to calculate the diagrams 
shown in Fig. [H 



Fv{t) 



QAl 



(t-4m2)J(t) 
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To calculate the matrix element ([H]) in a finite volume one has to evaluate the same diagrams, 
Fig. [H but with the propagator ([6]). We notice that the field renormalization is absent both 
in a infinite and finite volume. Calculating the tadpole diagram in Fig. [1] (a) one notices 
that the term with n = in the sum in Eq. ([6]) gives rise to the infinite volume result, while 
the remaining terms give corrections which vanishes at L ^ oo. It is convenient to carry out 
a similar rearrangement for the diagram (b) in Fig. [T] also. Schematically it amounts to the 
following transformations [22] 

j d'^xe'P''G'ix-y) = J2 j d'^xe'^^'Gooix - y + nL)G{x - y) = j d'^xe"'''Goo{x - y)G{x - y) 

= J d'^xe'P'^G^ix - y)Goo{x -y) + J^ J d^xe'^^Gooix - y)G^{x -y + nL). (13) 

ft 
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Again the first term in the last hne of Eq. (fT3l) represents the infinite volume contribution 
while the last one gives 1/L corrections. As usually the sign means that the term with 
n = in the sum has to be omitted. The integration in last line of Eq. (fT3!) goes over infinite 
space and involves infinite volume pion propagators ([7]). We notice also that each term in the 
sum (fT5]) is UV finite. One also has to take into account that due to the imposed boundary 
condition all spatial momenta are quantized, q = {2iT/L)n, where n = (ni,n2,n3) is a vector 
with integer components. Then, after lengthy but straightforward calculations one obtains the 
following answer for 1/L corrections to the matrix element ([8]) 

{n\p')\0^'^'\'K''{p))i/L = ^le^^^i^P^" {A~B)~ P^'C + A'^(P5)) . (14) 
Here P = (0, P) and the formfactors A, B, C, are given by the following expressions 

^=ZX^5Z'/ dr]M^{A,r])Ki{m{r]),n), (15a) 



B 
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^'Ki(m,n), (15b) 



'^=^5Z'^'/ t^^M„(A,77)K2(m(77),r?), (15c) 



dririMl{A,ri)Ki{m{ri),n) , (15d) 



^ 2tt ^ 

where rn?{vi) = rn? — ^{1 — rf )A^ , = (Aq, A) and we put A = — k. The function M„, 

L 

and K are defined as follows 



M„(A, v) = W {-ir-^- cos(7rr7 n^k^) , MI{A, r/) = ^-^^!pJ^M^{A, r/) (16) 

and 

' Kim\n\V) , (17) 
\n\ J 

where Ki{x) is the modified Bessel function of the second kind. The sums in (JT5l) converge 
rapidly since Ki(x) ~ x~^^'^e~^ for large x. Since the Lorentz symmetry does not hold for 1/L 
corrections the matrix element (fl^ is parameterized by three formfactors instead of one in an 
infinite volume limit. Let us remark that the difference A — B vanishes in the limit A — > 
which ensures that the temporal component, = 0, of matrix element ([8]) does not receive 
corrections in this limit. 



tors, A^2fl^ and 



The first moment of the isoscalar pion GPD, "(a^), is parameterized by two formfac- 

dxH'='{x, e, t) = Aifit) + AeAifit) . (18) 
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They enter the matrix elements of the isoscalar operator, 
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S q'j^iD^q- -g^^qWq 



{^\p')mz)\ir'^{p)) = 5'^'[Ai='{P . zf + (A ■ zf) 



(19) 
(20) 



where z is auxihary hght-hke vector, D,, 



Df_,- Df,), S Qf,u = UQ/iu + Qufi) and we set 



Q{z) = z^ZyQ^" . We also notice, that throughout the paper the scalar product (a ■ h) implies 
the Minkowsky scalar product. 

In ChPT the operator (IT^ reads 



(21) 



In an infinite volume the formfactor A20 does not receive nonanalytic corrections at order 
0{q^), while the result for the formfactor A22 reads [2^ 



^22 °(^) — ^22^ ( 1 



(0) 



log 



m" 4 t + 2m 



2 + 3 



J{t) 



(22) 



We recall also that = —A20 /A = 0.20- 

The volume corrections to the matrix element fl20|) comes from diagram (b) in Fig. ([1]) only, 
the contribution due to the tadpole diagram. Fig. [I](a), vanishes. The calculation gives the 
following result for the matrix element 



{Ap')\Q2i^)\^''{p))l/L =52,0 [t-^] T{Z) 



m 



Tiz] 



where T{z) = z^Zi,T^^ and 



D 



-y' t 



4A2 — J 1 



{A-zyD-{A-z){zS) + -z'C) 



dri (1 - ri^) M„(A, r/) Ko(m(r/), n) . 



(23) 
(24) 



(25) 



It is interesting to notice that both the infinite and finite volume corrections vanish at 
2t = . To make the structure of corrections more transparent we write down the tensor T^'^ 
in components 



T 



00 



{Ay-^A'' ]D + ^{AS) + ^C, 



1 



-lOi rpiO 



A° ( A'D 



=A^A^D + - [A'S'' + A^S') + ^ - {AS) 



(26a) 
(26b) 
(26c) 
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Thus parametrization of the volume corrections for an isoscalar operator at one-loop level 
require only one additional function, D. 

5. The lowest moments of chiral odd GPDs, E^^, are related to the matrix elements of 
the operators Q\ = go"AuT"g and Qxlo = A S qaxuiDpq [19]. The symbols S and A denote 

(anti)symmetrization in the corresponding indices. The pion matrix elements of these operators 
are parameterized as follows 

{^'(pmyip)) =^6«^^i?A,i??i,o(t) , (27) 

Wm^ijAp)) =^'^A SR,,P,B-,^,it) , (28) 

^'^ A/1 pp 

where Rx^ = (-PaA^ — AxP^)/m.„. The formfactors -Bf „o equal to the moments of the 
chirally odd GPD E^^, ^-BJ,no(^) ~ jlidxx^~^ET^n{x,C,,t)^. In ChPT these operators take 
the form (we follow the notations of Ref. [20] ) 

Qlp =\F%Tifl tr {u^t'^u^ + ut'^u) [u\ u^] , (29) 
^V.P =\F'bT2,, tr {u' + («t)2 j [u\2zW^)u'^ - (A - /.)] , (30) 

where = U, 

= i {u^d^u - ud^u^) , Vf,X = d^X + [F^, X] , = ^ {u^dpU + ud^u^) . 

In Eq. (!30l) it is implied that one has to subtract all traces and implement the indicated 
(anti)symmetrization 

2^V:P _^ 1 ^2T^^''^ + T^^''^ — T^'^''^ — g^^T^'^^ + g^'^T^'^^^ (31) 



The formfactors B^^^^it) at order O(g^) (in an infinite volume) were calculated in [20] 

Bhoit) =B?ll (l + ^ log ^ + Mt, m')^ + • • • , (32) 
^?,io(^)=5S(l-^log^)+..., (33) 

where the dots stand for analytic terms and -B^io = "^tt&ti.o, 

Calculating the diagrams in Fig. [1] with the insertion of the operators and Q^^^ ^ one 
obtains the following result for the volume corrections to the matrix elements 



{^"ip')\Qlp\^\p))i/L =^e"'' ^1 {P'^'' - ^^P^) ( ^ - Y ) " ^^'^"^ " ^'^'^ ^ \ ■ ^^^^ 
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We recall that P° = and P* = P*. Let us notice that 1/L correction to the matrix element 
of the chirally odd isovector operator (134|1 involves the functions A, B and C which enter the 
volume correction to the chirally even isovector operator (fT^ . 

In its turn the volume correction to the chirally odd isoscalar operator 



2 

m 



E-[t- — \C 



(35) 



involves one new function E, 

E=^Y^n^^,{m,n). (36) 

^ n 

We also introduced the notations = (0, S), where S is defined by Eq. ( ]15dp . Of course, the 
matrix element (135!) has to by symmetrized according to (13T|) . 

6. In this Letter we have presented the results of the calculation of volume corrections to 
the matrix elements of the pion operators of lowest spin. We found that the 1/L corrections 
to the four operators we considered involve (due to Lorentz symmetry breaking) 11 different 
kinematical structures with coefficients parameterized by six functions {A, B, C, D, E, S^). Our 
results should be helpful to improve the procedure of the extracting of GPD moments from 
lattice data. 
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